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ABSTRACT 

Using hydro dynamical simulations, we show for the first time that an episode of star 
formation in the centre of the Milky Way, with a star-formation-rate (SFR) nsj Q ^ ^ Mo 
yr ^ for ^ 30 Myr, can produce bubbles that resemble the Fermi Bubbles (FBs), when 
viewed from the solar position. The morphology, extent and multi-wavelength observa¬ 
tions of FBs, especially X-rays, constrain various physical parameters such as SFR, age, 
and the circumgalactic medium (CGM) density. We show that the interaction of the 
CGM with the Galactic wind driven by star formation in the central region can explain 
the observed surface brightness and morphological features of X-rays associated with 
the Fermi Bubbles. Furthermore, assuming that cosmic ray electrons are accelerated in 
situ by shocks and/or turbulence, the brightness and morphology of gamma-ray emis¬ 
sion and the microwave haze can be explained. The kinematics of the cold and warm 
clumps in our model also matches with recent observations of absorption lines through 
the bubbles. 

Key words: Galaxy: center - Galaxy: halo - ISM : jets and outflows - galaxies: star 
formation 


1 INTRODUCTION 

Galactic outflows are most easily seen in different wave¬ 
bands in starburst galaxies such as M82 and NGC253 (e.g., 
McCarthy 1987). However, the discovery of the gamma-ray 
emitting Fermi Bubbles (FBs) and the associated multi¬ 
wavelength features in the center of Milky Way (Su et al. 
2010) suggest that galactic outflows may be quite common, 
occurring in Milky-Way and lower mass galaxies frequently 
over their life-times. Supernova feedback is invoked to globally 
suppress star formation in halos less massive than 10^^M© 
(Shankar et al. 2006; Sharma & Nath 2013). However, it is 
well-known that isolated supernovae fizzle out in less than a 
Myr, requiring overlapping supernovae to form a superbubble 
(e.g., Sharma et al. 2014). Thus, galactic outflows driven by 
overlapping supernovae are likely to be ubiquitous. 

In this paper we present a hydrodynamical model for ex¬ 
plaining the multi-wavelength (gamma-rays. X-rays and ra¬ 
dio) morphology and brightness of FBs. We assume a super- 
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novae (SNe) driven Galactic wind as the driving mechanism 
for FBs. 

Several observations point toward the existence of a 
gaseous outflow from the centre of Milky Way. An enhance¬ 
ment in the diffuse soft X-ray emission in the longitude range 
—20° ^ I ^ 35° with an emission scale height (in the south¬ 
ern hemisphere) of b ~ —17° suggests a large-scale flow of 
gas out of the disc (Snowden et al. 1995; Everett et al. 2008). 
This emission was modelled by Snowden et al. (1995) with a 
mid plane gas density Ue ~ 3.5x 10“^ cm“^ and temperature 
T ~ 4x 10® K. Observations by Almy et al. (2000) proved that 
at least half of the central emission comes from more than 2 
kpc from the Sun, and most likely lies near the Galactic cen¬ 
tre (see also. Park et al. (1997); Yao & Wang (2007)). Almy et 
al. (2000) took into account other components (stellar, extra- 
galactic), and improved the model density and temperature 
to Ue ~ 10~^ cm“® and T ~ 8.2 x 10® K. Interestingly, this 
emission was predicted from a model of cosmic ray driven 
Galactic outflow by Breitschwerdt & Schmutzler (1994). In 
fact, using mid-infrared (8.3/im) and ROSAT (1.5keV) ob¬ 
servations, Bland-Hawthorn & Cohen (2003) first showed the 
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existence of a biconical Galactic outflow. They also speculated 
about the existence of projected x-ray bubbles on the both 
sides of the galactic plane extending up to ~ 80° in latitude. 

The discovery of 7 -ray bubbles in the similar part of the 
sky, known as the Fermi Bubble, has given spurt to exploring 
the high energetic implications of a Galactic outflow. These 
twin bubbles, extending up to ~ 50° (~ 8 kpc in height) 
above and below the Galactic centre, are marked by 7 -ray 
emission with a remarkably uniform surface brightness and a 
(dN/dE ~ E~‘^) spectrum that is harder than the emission 
from the disc (Su et al. 2010). 

The X-ray and 7 -ray features also coincide with emission 
features in other wavelengths, such as the microwave haze 
found by WMAP and Planck (Finkbeiner et al. 2004; Planck 
Gollaboration 2013) and the polarized radio lobes seen at 2.3 
GHz (Garretti et al. 2013). Incidentally, Lockman (1984) had 
noted a HI hole in the inner Galaxy. These morphological 
similarities, to the extent of the edges of the features in dif¬ 
ferent wavelength almost coinciding with each other, suggest 
a common physical origin. 

Several models have been proposed to explain the FBs. 
As far as 7 -ray emission mechanism is concerned, there re¬ 
mains an uncertainty whether the inverse Gompton scatter¬ 
ing of cosmic microwave background photons by relativistic 
electrons is the source (Su et al. 2010) or the interactions of 
high energy protons with protons in the medium (Grocker & 
Aharonian 2011). The high energy electrons or protons can 
either be accelerated in situ by internal shocks and turbu¬ 
lence (Mertsch & Sarkar 2011), or advected from the disc. 
Outflows triggered by star formation in the Galactic centre 
(GG) region (Grocker 2012 ; Lacki 2014) and by the black hole 
at the GG (Guo & Mathews 2012 ; Yang et al. 2012 ; Mon et 
al. 2014) have been proposed for the dynamical origin of the 
FBs. 

The AGN-based models (both jet and wind driven) gen¬ 
erally consider a shorter age (< few Myr) for the FBs because 
the inverse-Gompton cooling time (due to up-scattering of 
starlight) for ~ 100 GeV electrons responsible for gamma ray 
emission is a few Myr (Guo & Mathews 2012; Yang et al. 
2012). The speed required to reach ~ 10 kpc in 1 Myr is 
~ 10 ^ km s“^, achievable by relativistic jets slowed down by 
the hot circumgalactic medium (GGM). The power required 
for inflating young FBs is much higher and the outer shock is 
much stronger (Guo & Mathews 2012; Zubovas & Nayakshin 
2012), and X-ray emissivity and temperature much higher 
than what is observed (Kataoka et al. 2013). The SNe-driven 
models of FBs consider them to be long lived (> 10 Myr), 
and thus the injected power is smaller. In fact, Grocker et 
al. (2014b) suggest a hadronic origin for gamma ray emission 
and consider the FBs to be steady features older than few 100 
Myr. The outer shock is weaker in the SNe-driven models and 
the temperature and emissivity jumps are modest, consistent 
with X-ray observations. 

While the AGN jet and wind models have been explored 
numerically, simulations of a SNe driven model for FBs have 
not yet been carried out. The dynamical modelling is lim¬ 
ited to simple arguments invoking a steady wind, termination 


shock, thermal instabilities, etc. (Grocker et al. 2014b; Lacki 
2014). A realistic SNe-driven wind is expected to be affected 
by disc stratification and the presence of a GGM. Moreover, 
thermal and Rayleigh-Taylor instabilities are expected to mix 
the hot bubble gas with the halo gas. We capture these com¬ 
plex, time-dependent, multidimensional effects in the hydro- 
dynamic numerical simulations presented in this paper. 

In our model, the FB is a time dependent phenomenon 
and is currently expanding. Our goal is to study the time 
dependent signatures of a star formation triggered Galactic 
wind, and to identify various features observed in different 
wavebands ( 7 -rays, X-rays, microwave and radio) with var¬ 
ious structural features of a Galactic wind. In doing so, we 
pay particular attention to projection effects from the vantage 
point of the solar system in the Galactic disc. 

The paper is organised as follows. In section 2, we discuss 
the initial set up, the simulation settings and the parameters. 
The morphology and the importance of projection effects in 
the observations of the FBs are described in section 3. In sec¬ 
tion 4, we describe the X-ray, microwave and 7 -ray emission 
from the FB. Section 5 points out some kinematics aspects 
of the cold/warm clumps. We discuss the implications and 
improvements of our work in section 6 . Finally, in section 7, 
we summarise the main conclusions of this paper. 


2 SIMULATION 
2.1 Initial set up 

The details of the initial set up for milky way (MW) type 
galaxies are given in a previous paper Sarkar et al. (2015). 
However, for the sake of completeness, we briefly discuss the 
set up below. 

In our set up, we consider two gas components, a warm 
component (T = 4 x 10^ K, including the contribution from 
non-thermal pressure) representing the disc gas, and a hot 
component (T = 2.5 x 10® K) representing the extended 
circum-galactic medium (GGM). Since the warm gas repre¬ 
sents the disc, we also consider azimuthal rotation for this 
component. The hot GGM, however, is considered to be non¬ 
rotating. 

These two gas components are considered to be in steady 
state equilibrium with background gravitational potential of 
the stellar disc and dark matter (DM). For the disc, we use 
the Miyamoto & Nagai potential (Miyamoto & Nagai 1975) 

^aUR,z) = - {a,b^0) (1) 

+ (a + v/i 2 -:p ^)2 

where a and b are the model parameters representing the scale 
length and the scale height of a disc of mass Mdisc respectively, 
and, R and z are the cylindrical coordinates. For the dark 
matter, we use a modified form of NFW profile (Navarro et 
al. 1996) introducing a core at the center. The modified form 
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parameters 

values 

Mvir(M©) 

IO 12 

■Wdisc(M0) 

5 X lO^o 

Thalo (K) 

2.5 X 106 

rvir (kpc) 

258 

c 

12 

a (kpc) 

4.0 

b (kpc) 

0.4 

d (kpc) 

6.0 

^disc (-^o) 

1.0 

^halo {^q) 

0.1 

Pd(0, 0) (mpCm-3) 

3.0 

P/i(0, 0) (mpCm-3) 

2.2 X 10-3 


Table 1. Parameters used in our simulations. Hot gas central den¬ 
sity ph (0, 0) is obtained after normalising the total baryonic content 
(stellar plus gaseous) to 0.16 of Mvir, consistent with the cosmic 
baryonic fraction. While exploring the parameter space, we make 
this assumption flexible. 


of the potential is given as 


^dm(^, z ) 


( GM^ir '\ log(l + VR^ + z^+dyrs) . 

V/(c)rJ ^ 


( 2 ) 

where /(c) = log(l + c) — c/(l + c) with c = r^iv/rs as the 
concentration parameter and d is the radius of the core which 
gives a finite DM density at the centre. Cvir and Fs are, re¬ 
spectively, the virial radius and scale radius for a DM halo of 
mass Mvir- 


The steady state density distribution in a combined po¬ 
tential $(i?,z) = $disc(^, + ^DM(i?, ^) for the warm gas 
can be written as 


pd{R,z) = pd(0,0)exp(—h 2 ;) _ $(o, 0) 

-f($(fl,0)-$(0,0))]) , (3) 

and for the hot CGM, 

=/5h(0,0) exp (—- $(0,0)]) , (4) 

V ^sh / 


where, Pd(0, 0) and ph{0, 0) are the warm and hot gas central 
densities and Csd and Csh are the isothermal sound speeds of 
the warm disc and the hot CGM, respectively. Here, / is the 
ratio of the disc gas rotation velocity and the stellar rotation 
velocity at any R and taken to be a constant (= 0.95). The 
density of a given location is, therefore, pd-\- ph- A full list of 
model parameters is given in Table 1. 


2.2 Code settings 


and from ^ = 0to^ = 7r/2in theta direction using uniform 
grids. This implies that the disc lies on the Q — 'k j2 plane and 
our simulation box includes the first quadrant of the 2D slice 
taken along the d plane of our Galaxy. 

In our simulation, we express the temperature as T ^ pip 
which includes the hot gas pressure in addition to the 4 x 10^ 
K gas pressure inside the disc. The effective temperature of 
the disc is large enough to induce strong cooling unlike the 
warm gas at T = 10^ K gas. In reality the disc gas is always 
being heated by the supernovae and other processes. Since 
we are interested in the Galactic wind and not the disc ISM, 
we constrain the cooling of the disc material (but not the 
injected material) to be zero for a height less than 1.2 kpc 
above the disc plane. A more detailed description about the 
code implementation can be found in Sarkar et al. (2015). 


2.3 Injection parameters 


The mechanical luminosity of a starburst activity can be writ¬ 
ten as 


in40 -1 

10 ergs eo .3 


/ SFR \ 

VO.lMoyr-iJ ’ 


( 5 ) 


where, eo .3 is the thermalisation efficiency in units of 0.3 
and SFR is the star formation rate. Here we have consid¬ 
ered Kroupa/Ghabrier mass function, for which there is ~ 1 
SN for every 100 M© of stars formed. 

As we show later, the morphology and X-ray emission 
properties of FBs depend mostly on the combination of C 
and the GGM gas density. After scanning through various 
combinations of these two parameters, we show later (in §4.1, 
Figure 4) that a fiducial combination of £ = 5 x 10‘‘° erg s“^ 
and pho = 2.2 x 10 “® cm“® best matches the observations. 
The implied star formation rate, according to eqn 5, is ~ 
0.5 M© yr“^ (considering eo .3 = !)• The current rate of star 
formation in the central molecular zone of Milky Way is of 
order 0.1 M© yr“^. Mid-infrared observations by Yusuf-Zadeh 
et al. (2009) have led to an estimate of SFR ranging between 
0.007-0.14 M© yr“^, over the last 10 Gyr. Observations of 
young stellar objects in the central molecular zone (GMZ) in 
the 5-38 pm. band with Spitzer allowed Immer et al. (2011) 
to estimate a SFR of ~ 0.08 M© yr“^. The diffuse hard X- 
ray emission in the Galactic centre region was used by Muno 
et al. (2004) to estimate an energy input of ~ 10^° erg s“^. 
However, the star formation activity in the central region of 
the Galaxy is likely to be episodic. Our fiducial SFR, averaged 
over the last several tens of Myr, is therefore not unreasonable 
although it is a few times larger than the current SFR. 

The mass injection rate has been taken as (Leitherer et 
al. 1999) 


Minj = 0.3 SFR. 


( 6 ) 


We use the publicly available code PLUTO (Mignone et al. 
2007) for our hydrodynamic simulations. We perform the sim¬ 
ulations in 2D spherical coordinates assuming axi-symmetry 
around ^ = 0. The simulation box extends from rmin = 20 pc 
to rmax = 15 kpc in radial direction using logarithmic grids 


In our fiducial simulation, the considered mechanical lumi¬ 
nosity, £ = 5 X 10^° erg s“^, corresponds to SFR = 0.5 M© 
yr“^ and therefore Minj = 0.15 M© yr“^. 

We inject this mass and energy in density and energy 
equations inside a region of r ^ rinj (60 pc). The injection 
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Name 

C (erg 

:s-i) 

PhQ 

(> 

mp 

cm ^) 

SI 

1.0 

X 

1040 

0.5 

X 

10- 

-3 

S2 

1.0 

X 

1040 

1.0 

X 

10- 

-3 

S3 

1.0 

X 

1040 

3.0 

X 

10- 

-3 

S4 

2.0 

X 

1040 

0.5 

X 

10- 

-3 

S5 

2.0 

X 

1040 

1.0 

X 

10- 

-3 

S6 

2.0 

X 

o 

o 

1—1 

2.0 

X 

10- 

-3 

S7 

2.0 

X 

1040 

3.0 

X 

10- 

-3 

S8 

4.0 

X 

o 

o 

1—1 

0.7 

X 

10- 

-3 

S9 

5.0 

X 

1040 

1.1 

X 

10- 

-3 

SIO* 

5.0 

X 

o 

o 

1—1 

2.2 

X 

10- 

-3 

Sll 

6.0 

X 

1040 

3.0 

X 

10- 

-3 

S12 

1.0 

X 

o 

1—1 

1.1 

X 

10- 

-3 

S13 

1.0 

X 

10"^! 

2.2 

X 

10- 

-3 


Table 2. The list of runs showing the injected mechanical lumi¬ 
nosity and the central density in column 2 and 3, respectively The 
fiducial run (SIO) has been pointed out by an asterisk in the list. 


10000 100000 1e+06 1e+07 0.0001 0.001 0.01 0.1 



14 12 10 8 6 4 2 0 2 4 6 8 10 12 14 


Figure 1. Snapshot of density (right panel) and temperature (left 
panel) contours at 27 Myr for our fiducial run (SIO). The wind 
structure has been pointed out by different labels, from outside to 
inside as, CGM: circumgalactic medium, FS: forward shock, CD: 
contact discontinuity, SW: shocked wind and FW: free wind. 


rates can then be written as 


and 


2 C 

3 (47r/3)rf„j 

Minj 

(47r/3)rf„j 


( 7 ) 

(8) 


where, p is the pressure. A full list of all the runs is given in 
Table 2. 


3 RESULTS: WIND BUBBLE MORPHOLOGY 

The result of an episodic explosive event at the centre of Milky 
Way would depend mainly on the rate of energy and mass in- 


IqIS .|q19 .|q20 .|q21 


Column density (cm'^) 



90 80 70 60 50 40 30 20 10 0 -10 -20 -30 -40 -50 -60 -70 -80 -90 


I (deg) 

Figure 2. Snapshots of column density from edge-on position but 
without projection effects (left panel) and Solar vantage point with 
projection effects (right panel), for the same physical parameters 
as in Figure 1. The boundary of our simulation box (15 kpc) cor¬ 
responds to an angle ~ 60° from a distance of 8.5 kpc, and shows 
up in the left panel. 


put (and therefore on the SFR), the distribution of density 
through which the bubble ploughs its way (the disc and CGM 
gas density profile) and the epoch under consideration. We 
fix these parameters based on the morphology of the result¬ 
ing bubble, in light of the observed morphology of the FBs, 
and the emission properties. Therefore, we first discuss the 
morphology. 

Figure 1 shows the colour-coded contours of density and 
temperature for our fiducial run, £ = 5 x 10^° erg s ^ at t = 
27 Myr (corresponding to 10^ supernovae over this time). 
The snapshot clearly shows the structure of a standard stellar 
wind scenario (Weaver et al. 1977). There is an outer shock (at 
a vertical distance of ~ 10 kpc), an enhancement of density 
in the shocked CGM/ISM and shocked wind region, near the 
contact discontinuity (at a vertical distance of 6-8 kpc), as 
well as the inner free wind region (below a vertical distance 
of ~ 6 kpc). The figure also shows a second reverse shock 
at height of ~ 2 kpc which arises because of the presence of 
two component density structure related to the GGM and the 
disc. 

Since we are at a distance of 8.5 kpc from the centre of 
the Galaxy, and the wind-cone extends ~ 4 kpc at a height 
of 5 — 6 kpc, much of the observed structure is influenced 
by geometrical projection effects. Figure 2 illustrates the idea 
by showing the map of column density as viewed from an 
edge-on vantage point from infinity, as well as its appear¬ 
ance from the point of view of the solar system. In order for 
the column density not to be dominated by the disc mate¬ 
rial, we have considered only the gas for which the total non- 
azimuthal speed is larger than 20 km s“^. From 

the edge-on position, the Galactic coordinates are computed 
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as I = tan“^(_R/8.5kpc), b = tan“^(z/8.5kpc)/ whereas, 
from Solar view point (right panel), we have considered the 
projection effects accurately (see Appendix A for details). In 
projection from the Solar system position, the bubble appears 
bigger in angular size. Note that we have used the axisymme- 
try property of our 2D simulations to get the projected maps 
presented in this paper. 

The difference between the left and right panels of Figure 
2 highlights the importance of taking projection effects into 
account when comparing the morphology of the simulated 
bubble with the observed FBs. With the projected column 
density map at hand, we can discuss the logic behind fixing 
the epoch of the phenomenon at 27 Myr. 

As explained below, the X-ray emission expected from 
the outer shock (shocked circumgalactic medium, CGM) is 
likely associated with the observed Loop-I feature in X-rays. 
This feature is also observed in soft 7 -rays. The location of the 
outer shock depends strongly on the time elapsed, and helps 
us to fix the time at 27 Myr. The radius of the outer shock 
in a constant luminosity-driven wind, according to Weaver et 
al. (1977), is given by IZ ^ 


10 kpc 


( ^ 

O.OOlmp 

t 

1 5 X 10^°erg s~^ 

P 

27Myr 


1/5 


( 9 ) 


matching the outer shock location in Figure 1. Moreover, with 
this choice, we find that the location of the contact discon¬ 
tinuity matches the edge of the FBs. This indicates that the 
emission in different bands coming from the FBs is created 
within the the contact discontinuity. In addition, as shown 
below, the morphology of emission in different wavebands re¬ 
markably matches the predictions based on this choice of time 
elapsed (namely, 27 Myr) and therefore, in turn, supports the 
idea that some part of the Loop-I feature is likely associated 
with the FBs. A point to note in Eq. 9 is that the outer shock 
radius depends more sensitively on time rather than SFR or 
the CGM density. 

While Eq. 9 is strictly valid only for an homogeneous 
and isotropic medium, and with isotropic energy injection, 
we expect it to be roughly valid, even with anisotropic AGN 
jets. Most AGN-based models consider a shorter age (~ 1 
Myr), which comes at the expense of a much larger mechanical 
power (up to 10^^ erg s“^, Guo & Mathews 2012; Yang et al. 
2012). The velocity of the outer shock is given as V ~ 371 jbt 


^ 200 km s ^ 


( ^ 

O.OOlmp 

lOkpc 

1 5 X lO^^erg s~^ 

P 



( 10 ) 

comparable to the sound speed in the hot CGM (~ 180 km 
s“^), implying a weak shock in case of £ = 5 x 10^° erg 
s“^ as seen in Eigure 1 . A more powerful AGN jet acting 
for 1 Myr with C ~ 10^^ erg s~^ will result in a very strong 
shock, ruled out by X-ray observations that show only a slight 


^ These formulae are valid only for R, 2 ; ^ 8.5 kpc , or equivalently 
Z,6<45°. 


enhancement of temperature and density across the EB edge 
as observed by Kataoka et al. (2013). 

Though we assume that the injection region is spheri¬ 
cal symmetric, a departure from this assumption does not 
change the qualitative/quantitative picture much. The effect 
of different injection geometries has been discussed in section 
6.4. 


4 RESULTS: EMISSION IN DIFFERENT 
WAVEBANDS 

We discuss the results of our calculation for the emission in 
different bands in this section, and compare with the observed 
features. Various emission mechanisms have been discussed in 
the literature for different bands - gamma-rays. X-rays, mi¬ 
crowave and radio, and most of the debate so far has centred 
around the 7 -ray radiation mechanism (hadronic or leptonic), 
whether or not particles are being advected from the disc or 
accelerated in situ. However, among the emission in differ¬ 
ent bands, the X-ray emission from thermal gas suffers the 
least from any assumptions regarding accelerated particles 
and magnetic fields. We, therefore, discuss the X-ray emis¬ 
sion first. 

Since we have estimated the age of the Fermi bubbles to 
be 27 Myr (as discussed in the previous section), we perform 
detailed analysis for the fiducial run (SIO; see Table 2 ), at 
t — 27 Myr in this as well as in all the following sections. 


4.1 X-ray 

Observationally, two limb-brightened X-ray arcs, called 
‘northern arcs’, are seen in the north-east quadrant adjacent 
to the FB. In the southern hemisphere, a ‘donut’ feature is 
observed. Then there is the Loop-I feature extending up to 
b ~ 80° and from 50° to —70° in longitude. The diffuse X-ray 
emission also shows a dip in intensity in the FB region (Su 
et al. 2010). Recently Kataoka et al. (2013) have scanned the 
FB edge to look for differences in the X-ray brightness. They 
found that the temperature does not vary across the edge but 
there is a 50% decrease of the emission measure (EM) when 
moving from outside to inside of the bubble. 

We show the surface brightness of X-ray emission from 
the simulated bubble in Figure 3, in the 0.7-2.0 keV band 
(ROSAT R6R7 band) considering the Mekal plasma model 
(from XPEG) for emission. While calculating the X-ray sur¬ 
face brightness, we also consider the contribution from an 
extended GGM where the hydrostatic state has been extrap¬ 
olated to 100 kpc. ^ 

From Fig 3, we find that (1) the diffuse emission has a 


^ For our isothermal CGM the density profile is centrally peaked 
so that surface brightness (SB) is dominated by the inner CGM, 
but halo gas extended out to 100 kpc makes a non-negligible con¬ 
tribution (50%) compared to when it is confined to the 15 kpc box. 
The contribution of extended halo depends somewhat on the CGM 
density profile at large radii which is observationally unconstrained. 
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Figure 3. Simulated X-ray emission map in 0.7-2.0 keV band for 
the fiducial run (SIO), over plotted with the observed edges of the 
Loop-I, northern arcs and the northern FB. The white circle rep¬ 
resents the region where we have compared the estimated emission 
measure with the observations mentioned in the text . 


dip at FB and extends to the Loop-I feature in the form 
of a parachute (this fixes the age of FBs in our model, 
as mentioned earlier), and roughly delineates the feature 
leaving aside the slight asymmetry. (2) The location of the 
two arcs roughly matches the enhanced brightness (between 
60° < b < 50°, at / ~ 0). (3) The surface brightness is 
6 X 10“^ erg s“^ cm“^ sr“^, also consistent with obser¬ 
vations Su et al. (2010) who found ~ 2 x 10“® erg cm~^ 
sr“^. Although the contrast of the X-ray dip in the figure 
shown here is a bit less than observed, we note that the final 
observed counts through an instrument will depend on the 
details of spectral modelling. It is clear from figure 1 that 
the gas inside the bubble has a temperature (~ 2 X 10® K), 
lower than the shell temperature (~ 3.5 X 10®). Therefore, 
the intensity when folded through an instrument to estimate 
counts will show a higher contrast. These simulated features 
reasonably match the observed structure in X-ray images. 
Since the intensity of the ‘X-ray parachute’ mainly depends 
on the background CGM density, we use the emission mea¬ 
sure (EM) of the parachute at (/, 6 ) (10,55)° to match the 

EM of N1 point as observed by Kataoka et al. (2013) (shown 
by the white circle in Eigure 3). Eigure 4 shows the estimated 
values of the emission measure (EM = f rddl) for 0.24-0.38 
keV gas compared to the observed value of 0.05 cm“® pc for 
different runs (see Table 2). The figure shows that only for 
the central densities of 2-3.5 x 10“® m^ cm“®, the estimated 
EM is close to the observed value. This, along with the sur¬ 
face brightness of the ‘X-ray parachute’ sets a constrain on 
the allowed background density of the CGM. 

The different straight lines in figure 4 represent differ¬ 
ent values of CIpho (this ratio determines the radius of the 
outer shock for a given time; see Eq. 9). These lines of con¬ 
stant Cjpho are also found to be crucial in determining the 
shape of the wind, and hence the projected shape within the 
contact discontinuity. Eor lower values of Cjpho, the opening 


Figure 4. Comparison of the estimated emission measure (EM) of 
the ’X-ray parachute’ with the observed value 0.05 cm“® pc. The 
filled circles represent the position of individual runs in parame¬ 
ter space (as mentioned in Table 2) and the colour of each point 
represents the value of EM/0.05. Different values of Cjpho in lO"^® 
erg s“^ nip ^ cm^ have been shown by different straight lines. The 
fiducial run has been shown by ’SIO’ in this Figure. 

angle of the wind is much smaller than observed in FB (we 
assume that the gamma-rays of FBs come from the free and 
the shocked wind, as we discuss later). For larger values of 
Pm ^ though the opening angle matches with the base of 
the FB, the extent of the wind (in 1) at high latitudes exceed 
the observed width of the bubble. However, the shapes aris¬ 
ing from the runs lying on CIpho = 2 x 10 ® x 10 ^® erg s“^ 
m” ^ cm® line have maximum similarity with the observed FB 
shape. 

Therefore, the constraint on pho (~2 — 3.5x10“®) from 
X-rays and the requirement for the FB shape leave us with a 
small parameter space in figure 4 which implies C ^ 5-7 x 10 ^® 
erg s“^. Since modelling of thermal X-rays is least uncertain 
as compared to the non-thermal radio and gamma-ray emis¬ 
sion, we consider C — 5 x 10^® erg s“^ and pho — 2.2 x 10“® 
m.p cm“® for our Galactic wind model parameters to calculate 
microwave and 7 -ray emission. 

4.2 Microwave Haze 

Microwave observations (23 GHz, with WMAP and Planck; 
Dobler & Finkbeiner 2008; Planck Gohaboration 2013) show 
emission from | 6 | < 35° region on either side of the plane, 
termed the ‘microwave haze’. Diffuse radio emission is also 
seen in the 408 MHz map (Haslam et al. 1982) where the 
emission traces the Loop-I feature. The 23-70 GHz emission 
spectrum shows a spectral index /3 = 2.56 (brightness temper¬ 
ature Th oc h'~^) which indicates the presence of an electron 
spectrum of spectral index x = 2.2 (Planck Gohaboration 
2013). The 2.3 GHz observation also reveals polarised lobes 
and ridges in both hemispheres. The polarisation level in the 
ridges is 25-31%. The ridges in the north-east quadrant co¬ 
incides with the FB edge and the x-ray shells, and it is found 
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that the magnetic held is aligned with the ridges (Carretti et 
al. 2013). The low-frequency emission extends westward be¬ 
yond the FBs in both hemispheres and the spectrum 2.3-23 
GHz spectrum becomes softer as we go away from the Galac¬ 
tic center. 

In order to estimate the emission from relativistic parti¬ 
cles, we assume that the particle (either hadrons or leptons) 
energy density is a fraction of the total energy density of the 
gas (internal or thermal energy as well as the energy density 
due to huid motion) . This is expected in the case of internal 
shocks and turbulence in the gas, and due to in situ acceler¬ 
ation of particles from these shocks. As Figure 1 shows, the 
strong shocks that are likely to accelerate particles are traced 
by the shocked wind material which is the region inside the 
contact discontinuity (GD). Therefore, in order to trace the 
freshly produced accelerated particles, we use a tracer in the 
simulation that tracks the injected wind material from where 
most of the microwave/gamma-rays are emitted. Also, in or¬ 
der to avoid the disc material along the line of sight, we dis¬ 
card from our analysis the gas with a non-azimuthal velocity 
less than 20 km s“^. 

We can estimate the microwave emission in our model 
assuming synchrotron emission and that the cosmic ray (GR) 
energy density is given by Ucv — Ccr , where Ugas = 
uth + Ukin is the total energy density of the gas as discussed 
above. The GR electron energy density is assumed to be 
Ucr,e — 0.05?Xcr as cxpectcd from the ratio (rne/mp)^^~^^^‘^ 
for X — 2.2 (see, e.g., Persic & Rephaeli (2014)). This hxes 
the electron spectrum, n{E)dE = HiE~^dE, where the nor¬ 
malisation constant k, is given by = UcT,e{x — 2) 

(assuming a lower cut-off of Lorentz factor ~ 1). The syn¬ 
chrotron emissivity, in the presence of a magnetic held B in 
the optically-thin limit, is (Eq. 18.18 in Longair (1981)). 


syn 

erg s“^cm“^Hz“^ 


X 


1.7 X 10“^^a(x)/^H^ 

1^6.26 X 10^® Hz^^ 


( 11 ) 


where, a(2.2) 0.1. For the magnetic held, we assume that 

the magnetic energy is also a fraction of the thermal energy 
and is given as ub = cb ^gas • 

We therefore have the volume emissivity per unit solid 
angle as 


jsyn 


erg s“^cm“^Hz ^sr 


_ o vy 1 n 0.8 1.8 

:j-- = 2.6 X 10 CcrCB P 


23 GHz 


( 12 ) 


where, we have taken Ugas = (3/2) p, and p = pth + 1/3 pv^ is 
the total pressure (thermal plus kinetic). 

After calculating the surface brightness at 23 GHz from 
the FBs with Ccr = cb = 1.0, we found it to be approximately 
15 times larger than the observed value of 800 Jy sr“^. This 
implies that 


0.8 

Ccr Cb 


1/15. 


(13) 
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Figure 5. 23 GHz synchrotron emission (surface brightness) map 
for C = b X lO"^® erg s“^ and = 2.2 x 10“^ m^ cm“^ with 
Ccr = 0.15 and es = 0.4. The upper colourbar shows the brightness 
temperature in mK and the lower colourbar shows the brightness 
in units of Jy sr“^. 



Note that we have an independent constrain on Ccr because 
these same particles will also emit 7 -rays. Assuming Ccr = 
0.15, we get a constrain on the magnetic energy density that 
CB — 0.4. This gives a magnetic held of strength B = 3-5 pG 
considering Ugas ~ 0.7-3.0x 10“^^ erg cm“^ inside the bubble; 
thus, 23 GHz emission comes from electrons with 'y fb 2 x 
10^. Notice that our estimate of magnetic held is somewhat 
lower than but consistent with other estimates (Su et al. 2010; 
Garretti et al. 2013; Grocker et al. 2014b). 

The surface brightness of the 23 GHz emission is shown 
in Figure 5 and is consistent with observations. We also notice 
that the emission hlls up the whole bubble volume which is 
consistent with recent observation. Planck has detected mi¬ 
crowave emission from the whole FB region, although the 
intensity is rather small above ~ 5 ^ 40° ( see hg 9 of Planck 
Gollaboration (2013)), consistent with our results, given the 
uncertainties. 


4.3 7 -ray 

Observations show two 7 -ray bubbles (1.0-50.0 GeV) on ei¬ 
ther side of the Galactic plane, being roughly symmetric 
about the plane. The northern bubble extends up to 6 < 50° 
and |/| < 25°, which is almost same for the southern bubble. 
Another limb-brightened 7 -ray feature extends up to 80° in b 
and ±70° in I in northern hemisphere and is known as Loop-I 
feature. The FB surface brightness is fairly uniform over the 
bubble and shows no limb brightening. The 7 -ray spectrum of 
the FB is also flat (dN/dE ~ E~^) and shows almost no soft¬ 
ening with increasing height. The Loop-I feature, however, 
has a softer spectrum, dN/dE ~ E~‘^'^, and this has led Su 
et al. (2010) to conclude that Loop-I is a part of the disc and 
has no connection with the bubble. 

In order to estimate the 7 -ray emission from the simu- 


© 0000 RAS, MNRAS 000, 000-000 














8 Sarkar, Nath and Sharma 


10 '® 10 '® 

y-ray emission (GeV s'^ cm'^ sr'^) 

90 - 

80 

70 

60 

S 50 
0 

?40 
30 
20 
10 
0 

90 80 70 60 50 40 30 20 10 0 -10 -20 -30 -40 -50 -60 -70 -80 -90 
I (deg) 

Figure 6. Hadronic 7 -ray emission map (surface brightness) as 
seen from the solar system location. Over plotted are the edges of 
the observed emission maps. 



lated bubble, we consider two possible emission mechanisms, 
hadronic and leptonic. Below we discuss them in detail. 


4.3.1 Hadronic emission 

In the hadronic model, cosmic ray (CR) protons undergo 
hadronic collisions with thermal gas protons and produce 7 - 
ray via pion decay. The volume emissivity for this emission 
in the Fermi-LAT band (IGeV-lOOGeV) can be written as 

3 1 

_ X _ X fhoi Ucr,p'fT'p CTpp KppC , 

where, Up is the gas proton number density, u cr,p ~ UcT f is 
the GR proton energy density, fhoi — 0.4 is the fraction of 
the total luminosity that is emitted in the Fermi-LAT band. 
The factor 3/2 corrects for the presence of heavy ions among 
the beam and target nuclei, app = 4 x 10“^® cm^ is the 
corresponding interaction cross-section and Hipp = 0.5 is the 
hadronic inelasticity (Grocker et al. 2014a). 

Figure 6 shows the emission map in 1-100 GeV from 
hadronic collisions assuming Ccr = 0.15 as discussed in the 
previous section. The surface brightness is fairly uniform and 
fills the observed region of FB. However, the average inten¬ 
sity is only < 1 % of the observed value ofl.4xl0“® GeV s“^ 
cm“^ sr“^ (Su et al. 2010) in this band. 

4.3.2 Leptonic 

In the leptonic model, GMB photons are inverse Gompton 
scattered by relativistic electrons and produce 7 rays.^ The 
required GR electron energy is 1-100 TeV to produce 7 rays 


^ We do not consider up-scattering of UV and IR photons because 
CMB dominates far away from the galactic disc. A proper mod¬ 
elling of the 7 -ray emission due to interstellar radiation field (ISRF) 
is beyond the scope of this paper. 


of 1-100 GeV. Therefore, the corresponding Lorentz factors of 
the electrons range from 7 = 2 x 10® to 2 x 10^. However, at 
such high 7 values the electron spectrum would likely suffer a 
cooling break because of synchrotron and IG losses. Assum¬ 
ing a typical age of the electrons to be tage = 1.5 Myr, and 
considering the magnetic field to be H = 4 /xG, as obtained 
in section 4.2, we get the Lorentz factor at the break to be 
76 = 10®. Typically in a simple steady-state model for the 
evolution of the relativistic electron distribution function, the 
cooling break occurs at a 7 for which the cooling time equals 
the age. However, here we are considering time-dependent 
particle acceleration in turbulence and internal/termination 
shocks, and the effective age of electrons can be much shorter 
than the FB age. 

Since 7 = 10® is close to the Lorentz factors needed for 
the leptonic emission to be in 1-100 GeV band, we consider 
a broken power law electron spectrum which has a spectral 
index xi — 2.2 (same as considered in synchrotron emission) 
below the break and the index drops by Ax = 1 after the 
break. The electron spectrum can be written as 


1 ( 7 ) = 




for 7( < 7 < 76 
for 76 < 7 < 7/1 


(15) 


where, X 2 = xi + 1 is the spectral index after the break, 7 ^ 
and are the lower and higher cut-off of the spectrum. The 
normalisation factor C can be written as 

n -1 


C = 


i^cr,e 

rUe C^ 




+ 


7 r^^ 


X 1 — 2 X 2 — 2 


xi — 2 


me c3 -)// 


(16) 


The output power per unit volume at an energy ei can be 
calculated as (Eq. 7.28a in Rybicki V Lightman (1979)) 


dE 3 ^ 

L-——= -cgtC e\ 
dV dtdei 4 


/ * (^) 


ri 

3 11 


7 - 1-7 


47^6 


djEjb 7 


ri 

d ih 


-X2-2 


47^6 


d'y 


(17) 


where, 

/ (x) = 2x log(x) + X H- 1 — 2x^, for 0 < X < 1 (18) 


and for seed photons with a blackbody spectrum at temper¬ 
ature Tcmb 5 


v{e) = 


Stt _£_ 

exp(€/fcsTcmb) - 1 


(19) 


Here we use xxcr,e = 0.05 x Ucr = 0.05ecr'?^^as with Ccr = 0.15 
(same as considered previously for synchrotron emission). The 
lower and the higher cut-off Lorentz factors are taken to be 
7 z = 1 and 7 /^ = 00 . Eqn. 17 can be numerically integrated to 
give the resulting spectrum, which is shown in Eigure 7. The 
spectrum shows a reasonable match with the spectra as ob¬ 
served by Su et al. (2010) and Ackermann et al. (2014). The 
figure also shows that the spectrum is consistent with the ob¬ 
servations for 76 ranging from 5 x 10® to 2 x 10® and therefore 
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Figure 7. Output spectra for leptonic 7 -ray emission (green 
dashed line). The blue empty squares and the magenta filled 
squares show the observed data points (Su et al. 2010 ; Ackermann 
et al. 2014), and the green (dashed) line shows the spectra calcu¬ 
lated by us for 75 = 10®. The plot also shows the spectrum for 
= 5x 10® (blue dotted line) and 7 ^ = 2 x 10® (red solid line) for 
comparison. Notice that we do not consider a high energy cutoff for 
the electron distribution, which can account for the lack of gamma 
ray emission beyond few 100 GeV. 
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Figure 8. Leptonic 7 -ray emission map at 10 GeV as seen from 
the solar system location. The red and black circles represent the 
regions for which the velocity histograms have been been shown in 
Figure 10. The white circle is the one where Fox et al. (2014) have 
UV absorption data. 


it is robust under small uncertainties in the magnetic field or 
age of the electrons. While the flux in 1-100 GeV decreases 
for a smaller 7 ?^, it can be boosted by the additional IC up- 
scattering of the ambient starlight. 

We also show the leptonic emission map at 10 GeV in 
Figure 8 . It shows a good match with the observed morphol¬ 
ogy of FBs. The surface brightness is also reasonably uniform 
over the region. Though the edge of the simulated bubble 
is not as smooth as observed, an introduction of magnetic 



01 23456789 10 

Distance in kpc 

Figure 9. Position-Velocity diagram of the gas parcels with T < 
2 X 10® K. The colourbar represents the density of the gas parcels. 

field in the simulation can potentially make the bubble edge 
smoother. 


5 RESULTS: KINEMATICS 

It is important to study the kinematics of FBs in order to 
infer their origin, as it can give us crucial information about 
the speed of gas inside and around the bubbles. Recently, Fox 
et al. (2014) have detected ultraviolet absorption features in 
cold (~ 5 X 10^ K) and warm gas (~ 10^ K) phases at line 
of sight velocities of —200, +130 and +250 km s“^ towards 
quasar PDS 456 (10.4°, 11 . 2 °). Using a simple model of bi- 
conical nuclear outflow, to obtain a line-of-sight velocity of 
~ —200 km s“^, they needed a cold/warm radial Galacto- 
centric outflow with velocity (ugsr) ^ 900 km s“^. This is 
essentially because of the radial outflow assumption and the 
low inclination of the quasar sight line. 

The velocity structure in our simulated FBs has a more 
complicated structure than the simple models studied by Fox 
et al. (2014). In our simulation, the cold/warm clouds are 
formed by thermal and Rayleigh-Taylor instabilities at the 
interface of the injected gas with the GGM (the contact dis¬ 
continuity). The clouds formed at the conical surface of the 
contact discontinuity sometimes fall back due to the gravity 
(essentially a fountain flow; Shapiro Field 1976). However, 
the clouds at the top of the cone keep moving away from 
the centre because of the wind ram pressure. The low lati¬ 
tude cold/warm gas can have a wide angle and a large (~ 100 
km s“^) line-of-sight velocity because the clouds are following 
non-radial trajectories (e.g., see the S 2 sequence of clouds in 
Fig. 12 of Sarkar et al. 2015). 

Figure 9 shows the posit ion-velocity diagram as seen from 
the Galactic centre. It shows that along with the positive 
velocity components of the warm gas, there are gas parcels 
which have negative velocities extending up to —100 km s“^. 
This infahing gas can contribute to the negative velocities as 
observed by Fox et al. (2014). 
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Figure 10. Velocity histograms of the gas along (0°,20°) in left panel and along (9°, 21.5°) in right panel. X-axis represents the line of 
sight velocity and the y axis represents the N(HII) corresponding to that velocity. The upper panel shows the velocity histogram for cold 
(T < 4 X 10"^ K) and for warm (4 x 10"^ < T < 2 x 10^ K), whereas, the lower panel shows for the hot gas (T > 10® K). The LOSs have 
been pointed out in Figure 8. 


In figure 10, we show the line-of-sight velocity (I’los) his¬ 
tograms of the cold, warm and hot gas along two different 
lines-of-sights (shown by the black and red circles in Figure 
8) that pass through the FB. We take into account the solar 
rotation velocity of = 220 km s“^ for this calculation 
(see Appendix A for more details). The central peak in all 
the histograms represent the stationary disc and halo gas. 
The upper panels clearly show that the ^'los for the cold and 
warm gas can reach up to —150 km s“^ and +200 km s“^. We 
also show the velocity-histograms of the hot (T > 10® K) gas 
in the lower panel of Figure 10. Though the hot gas velocities 
extend all the way from —200 to +600 km s^ for these two 
line of sights (LOS), the other LOSs show hot gas velocities 
extending from -500 to 700 km s~^. Notice that though the 
hot gas in our simulation has high velocity (~ 1000 km s“^) 
within the free wind region, the shape of the histograms dif¬ 


fer from that of the high velocity gas considered by Fox et al. 
(2014). This is essentially because of the non-radial flow of 
hot gas induced by Kelvin-Helmholtz instability. 

The above-mentioned results show that the kinematic 
signatures of our model are consistent with observations and 
one does not necessarily need cold/warm gas with velocities 
> 900 km 

6 DISCUSSION 
6.1 Diffusion of CRs 

In this paper, we have so far tracked CRs by using a tracer 
which confine them within the contact discontinuity (CD). 
In reality, CRs diffuse beyond the CD and produce extended 
7 -ray and synchrotron emissions. Since turbulence is not ex- 
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pected to be significant outside the CD, the in-situ acceler¬ 
ation of CRs becomes ineffective and the emission contains 
a signature of ageing in its spectrum. We can estimate this 
length scale over which the relativistic electrons diffuse before 
they lose their energy, by considering the diffusion coefficient 
to be (Gabici et ah 2007) 

For electrons of energy E' ~ 1 TeV moving in a magnetic 
held of R ~ 10 fiG (estimated from huid compression at the 
forward shock), the length-scale of diffusion in tage = 1.5 Myr 
(see §4.3.2) is cr = GDtage ~ 1.3 kpc. This implies that CR 
would diffuse up to ~ 10° beyond the CD at a height of 8 kpc. 
This extended emission can appear as 2.3 GHz radio lobe as 
observed by Carretti et al. (2013). 

6.2 Note on Loop-I 

The accelerated relativistic particles from the outer shock 
may also produce gamma-ray emission at the shock position. 
However, because of the absence of further acceleration mech¬ 
anisms {viz. turbulence) behind the shock, the particles will 
lose their energy and may have a spectrum that is different 
from that of FB. In case of leptonic emission at the outer 
shock, the electrons lose their energy after tage = 1.5 Myr 
giving rise to a faint gamma ray shell which can appear as a 
diffuse emission when viewed from the Solar vantage point. 
This may partially explain the observed emission from Loop-I 
feature. 

Incidentally, we note that recent observations by Acker- 
mann et al. (2014) (Fig. 13) have revealed a Southern counter¬ 
part of the Loop-I feature. This lends an additional support 
for the connection between the Loop-I and the FBs. 

6.3 Kinematics of the cold/warm gas 

The speed of the cold/warm gas is complex in nature as 
shown in Figure 9, ranging from —150 to +1000 km s“^. 
The density of the clouds, however, decreases with increas¬ 
ing velocity making them hard to detect. Moreover, very high 
velocity clouds may appear to be moving with low LOS ve¬ 
locity because of our vantage point. Therefore, even if star 
formation at the Galactic centre produces very high velocity 
clouds (VHVG), several factors can make them undetectable 
as VHVGs from the Solar vantage point. The +1000 km 
streak in the figure represents adiabatically expanding free 
wind. Because of its cone-like geometry, the velocity disper¬ 
sion along a line-of-sight can become large enough for a given 
column density that any absorption feature corresponding to 
the free wind may not be visible. 

6.4 Effects of the injection geometry and CGM 
rotation 

Although we have considered the injection region to be spher¬ 
ically symmetric with radius rinj = 60 pc, star formation can 


occur in a region of complicated geometry as observed by 
Molinari et al. (2012). To understand the effect of different 
injection geometries, we have carried out the following set of 
runs for our hducial value: 

• spherically symmetric, with rinj = 40 pc. 

• spherically symmetric, with rinj = 100 pc. 

• axisymmetric (about R = 0 axis) disc-like injection re¬ 
gion, with a radius R = 110 pc and midplane to edge height 
h = 42 pc (Lacki 2014). 

• axisymmetric ring-like injection injection region, with in¬ 
ner radius Rin = 70 pc, outer radius Rout = 100 pc and mid¬ 
plane to edge height /i = 50 pc. 

• spherically symmetric, with rinj = 60 pc, same as the 
hducial run but this time we have considered that the halo 
gas is also rotating with a speed equal to 10% of the stellar 
rotation speed at z = 0 at that R. 

The other parameters, namely ecr,es,7b(= 10®), have 
been kept hxed to those values mentioned in the text for the 
hducial run. The projected leptonic 7-ray emission maps at 
10 GeV for all the cases (at t — 27 Myr) have been shown in 
hgure A2. The hgures show that apart from a slight corruga¬ 
tion of the FB edge, the morphologies and intensities match 
quite well with each other. Therefore, the conclusions in this 
paper are not affected signihcantly by our assumptions re¬ 
garding the geometry of the star formation region and halo 
rotation. 

6.5 Effects of CR and magnetic pressure 

In our model, the magnetic and the GR energy densities are, 
respectively, 0.4 and 0.15 times the gas energy density. The 
addition of these forms of energy would surely increase the 
total energy content and hence the pressure inside the bub¬ 
ble. More so, because GRs and magnetic helds do not suffer 
signihcant radiation loses. Therefore, for the same bubble en¬ 
ergetics the required SNe energy injection, and consequently, 
the SFR will decrease approx, by a factor of 1.5 if we consider 
a strong coupling of these nonthermal pressures with the gas. 
This can bring down the required SFR to 0.3 — 0.4 M© yr“^. 


7 SUMMARY 

In this paper, we have presented the results of numerical sim¬ 
ulations of SNe driven outflows in our Galaxy, taking into ac¬ 
count a gaseous disc {T cb 10^ K) and halo gas (T = 2.5 x 10® 
K), in order to explain the origin of the Fermi Bubbles and 
multi-wavelength features related to it. We injected contin¬ 
uous SNe energy at the Galactic centre for 50 Myr. We as¬ 
sumed in situ acceleration of relativistic particles inside the 
contact discontinuity. Our model can explain the gamma-ray 
emission and microwave haze as coming from the interior of 
the contact discontinuity via leptonic and synchrotron emis¬ 
sion respectively. In addition. X-ray is emitted mostly by the 
shocked GGM. Given our vantage point at the Solar position, 
we considered the projection effects properly to calculate the 
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morphology of the bubble. In order to understand the dynam¬ 
ics of the bubble, we have also studied the cold, warm and 
hot gas kinematics. 

In our analysis, we have used the observed surface bright¬ 
ness in different bands in order to constrain the background 
CGM density, the star formation rate at the Galactic centre 
and the magnetic field inside the bubble, self consistently. We 
summarise the main results of our paper below. 

• The X-ray emission appears to have a parachute-like 
structure with a dip in intensity within the boundaries of 
the Fermi Bubbles. The surface brightness of the parachute is 
comparable to the observed value only if the central GGM gas 
density is taken to be 2-3.5 x m^ cm“^ and an extended 
GGM up to 100 kpc. Gonsidering the morphological aspects 
of the Fermi Bubbles along with the GGM, the injected me¬ 
chanical luminosity is found to be 5-7 x 10^° erg s ^ which 
corresponds a SFR 0.5-0.7 M© yrF 

• Assuming that a relativistic electron population of spec¬ 
tral index x — 2.2 gives rise to the microwave haze via syn¬ 
chrotron emission, we estimated the magnetic field inside the 
bubble to be 3-5 /xG which is little lower but consistent with 
other estimations. This electron population diffuses out from 
the contact discontinuity and produces polarised radio emis¬ 
sion as observed. 

• Gonsidering the above constrained GGM, SFR and the 
magnetic field, the 7-ray emission from the region inside con¬ 
tact discontinuity appears to have the shape and brightness 
comparable to observations. 

• The speed of the cold (T < 4 x 10^ K) and warm 
(4 X 10^ < T < 2 X 10^ K) clumps can vary from —150 km 

to +200 km s“^ (warm), whereas, the hot (T > 10® K) gas 
have a higher dispersion in their velocities which range from 
~ —500 to +700 km s“^. The kinematics of the cold/warm 
clumps appear to have the characteristics that are similar to 
recent observations. While Fox et al. (2014) argue for a large 
radial velocity outflow (e.g., associated with AGN-driven out¬ 
flows) because of a small angle between the radial direction 
and the line of sights through low latitudes. Figure 10 shows 
that we can obtain line-of-sight velocities consistent with ob¬ 
servations because of non-radial trajectories and the negative 
radial velocities of entrained cold clumps without requiring 
large radial velocity of cold/warm gas. 

To conclude, our work shows that star formation at the 
Galactic centre can give rise to the observed Fermi Bubbles 
and the multi-band structures related to it. Furthermore, 
modelling these structures can yield the basic Galactic pa¬ 
rameters such as the hot GGM density and magnetic field 
and opens up a new window to study high energy interac¬ 
tions in the Milky Way. 
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tion, can be written in terms of the Cartesian unit vectors as 

r = isin(0) cos(0) + j sin(0) sin(0) + kcos(0) 

6 = icos(^) cos((/)) + j cos(0) sin(0) — ksin(0) 

0 = -isin(0) + jcos(0). (A2) 

Therefore, the component of the actual velocity, ^ = rvr + Ovq + 
(pV(j), along the LOS is 

^los = . (A3) 

However, because of the solar rotation of on the plane, the 

actual LOS velocity along some (/, b) is 


'^los — 


^.L — V(j)^Q sin(/) cos(6) 

= — cos(b) cos{l) ^Vr sin(0) cos(0) + vq cos(6) cos(0) — sin(0)j 

— cos(b) sin(/) ^Vr sin(0) sin(0) + vq cos(6) sin(0) + cos(0)j 
+ sin(6) ^Vr cos{6) — vq sin(^)j —sin(Z) cos{b) . (A4) 


Figure Al. The geometry for calculating the projected maps and 
LOS velocity from the Solar vantage point. The coordinate system 
has its origin at the Galactic center (GC), whereas the Sun is situ¬ 
ated at a distance of Rq = 8.5 kpc from the GC along the X-axis. 
The Y-Z plane represents the plane perpendicular to the Galactic 
disc. The Galactic coordinates are marked as I and 6, the simula¬ 
tion coordinates are shown as r and 6 in the diagram and L is the 
distance of the point under consideration from the Sun. 
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APPENDIX A: PROJECTION EFFECTS 

In our analysis, we consider the Cartesian coordinates whose ori¬ 
gin is at the GC and the variables are calculated from the local 
standard of rest (LSR) and in terms of galactic coordinates (/, 5). 

Figure Al explains the working geometry. For a given (/,6), 
the unit vector along the line of sight can be given as 

L = — i cos(6) cos(l) — j cos{b) sin(Z) + k sin(6) . (Al) 
The unit vectors in spherical coordinates, used in the simula- 
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Figure A2. Effect of different injection geometry and CGM rotation for our fiducial run SIO (see Table 2) on the 7-ray emission at 10 
GeV. The geometries of the injection regions, from top-left to bottom-right, are i) spherically symmetric, with injection radius rinj = 60 
pc (the fiducial value), ii) spherically symmetric, with rinj = 40 pc, iii) spherically symmetric, with rinj = 100 pc, iv) axisymmetric, with 
disc-like injection zone, v) axisymmetric with a ring-like injection region and vi) spherically symmetric, with rinj = 60 pc and a rotating 
halo. The dot-dashed line shows the observed edge of the FB. 
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